Технические науки 


УДК 535.37 

ОПРЕДЕЛЕНИЕ ПОТЕНЦИАЛЬНОЙ ФУНКЦИИ МОЛЕКУЛЫ АбН 3 
НА ОСНОВЕ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ 

Ю.Б. Юхник, Е.С. Бехтерева*, Е.А. Синицын*, А. С. Булавенкова 

Томский политехнический университет 
Томский государственный университет 
Е-таіІ: уикЬпік@рНу$Т$и.ги 

Рассматривается задача определения внутримолекулярной потенциальной функции молекулы типа симметричного волчка на 
примере молекулы арсина АзН 3 . Для решения данной задачи разработан пакет программ на аналитическом языке МАРІЕ, по- 
зволяющий связывать между собой параметры потенциальной функции, заданной через симметричные, естественные (инвари- 
антные относительно изотопозамещения) и/или нормальные колебательные координаты. 


Известно, что внутренняя динамика молекулы 
в основном определяется внутримолекулярной по- 
тенциальной функцией [ 1 ] . Поэтому важным усло- 
вием понимания протекающих в молекуле процес- 
сов является получение корректной информации о 
последней. В связи со сказанным становится поня- 
тен интерес к количественному определению пара- 
метров потенциальной функции молекулы. 

Следует отметить, что такого рода проблемы ре- 
шаются в спектроскопии на протяжении уже нес- 
кольких десятков лет. Существуют два основных 
подхода ее решения: это использование аЬ іпіііо ме- 
тодов, однако точность таких расчетов всё ещё ос- 
таётся на 4-5 порядков хуже экспериментальных 
точностей в определении тех наблюдаемых вели- 
чин, которые могут быть получены на основе ин- 
формации о внутримолекулярной потенциальной 
функции. Другим подходом являются полуэмпири- 
ческие методы, которые основаны на использова- 
нии экспериментальной информации о молеку- 
лярных спектрах высокого разрешения инфракрас- 
ного и микроволнового диапазонов. В отличие от 
аЪ іпіііо методов, полуэмпирические методы позво- 
ляют получать информацию о параметрах потен- 
циальной поверхности с существенно более высо- 
кой точностью. 

Что касается молекулы арсина, и в частности её 
потенциальной функции, интерес к ней вызван 
несколькими причинами. С одной стороны, изуче- 
ние спектров арсина важно для прикладных задач 
астрофизики, поскольку они были обнаружены в 
атмосферах планет-гигантов Сатурна и Юпитера 
[2]. С другой стороны, такие объекты интересны 
чисто с теоретической точки зрения, т. к. все они 
легкие пирамидальные молекулы. Как следствие, в 
их колебательно-вращательных спектрах наиболее 
ярко проявляются все возможные эффекты и осо- 
бенности, которые могут присутствовать в спек- 
трах других молекул такого типа. Дополнительный 
интерес к исследованию спектров этих молекул 
вызван тем, что они удовлетворяют так называемо- 
му приближению локальных мод [3]. 

Заметим, что потенциальная функция АкН 3 
определялась лишь аЬ іпіііо методами [4] , либо от- 
дельные параметры (связанные с валентными ко- 


лебаниями) определялись на основе эксперимен- 
тальной информации в модели локальных мод [5]. 
Поэтому определение параметров потенциальной 
функции арсина на основе экспериментальных 
данных является актуальной задачей на сегодняш- 
ний день. 

В данной работе был проведен анализ центров 
колебательно-вращательных полос молекулы А«Н 3 
и на этой основе оценены квадратичные силовые 
параметры арсина, а также колебательные параме- 
тры ангармоничности и резонансные постоянные, 
позволяющие воспроизводить и предсказывать 
центры колебательно-вращательных полос с высо- 
кой точностью. Оценки были выполнены на осно- 
ве как имеющейся в литературе эксперименталь- 
ной информации [6-8] , так и на основе наших дан- 
ных, опубликованных в работе [9]. 

Хорошо известно [10], что все спектроскопиче- 
ские параметры молекулы являются функциями 
параметров равновесной конфигурации молекулы 
г‘ ш и констант форм колебаний І ш . Таким образом, 
становится понятным, что от точности их опреде- 
ления зависит точность всех дальнейших вычисле- 
ний. Поэтому на первом этапе, зная с высокой точ- 
ностью равновесную длину связи г‘ и равновесный 
угол а е , мы определили равновесные координаты 
ядер и константы форм колебаний молекулы А5Н 3 . 

Первые определялись из условия совпадения 
начала молекулярной системы координат с цен- 
тром масс ядер молекулы и условия совпадения 
осей молекулы с главными осями инерции 

Х'АѵТѵ а = 4 =°> а *Р> (1) 
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где N=1, 2, 3, 4 - число ядер в молекуле А$Н 3 , т к — 
масса соответствующего ядра, г‘ Ма - значение рав- 
новесных координат ядер. Отметим, что для упро- 
щения решения данных уравнений следует вос- 
пользоваться свойствами симметрии молекулы. 
Решение системы уравнений (1) дает 
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где г‘=\, 51101 А - значение равновесной длины 
связи, а'=92,0690° - равновесный угол между свя- 
зями [9] ,т и ,М- массы ядер водорода и арсина, со- 
ответственно. 

Зная равновесные координаты ядер, далее можно 
найти константы форм колебаний молекулы А$Н 3 . 
Заметим, что из анализа только симметричных 
свойств молекулы можно выразить все ЗА(3/Ѵ-6)=72 
параметров І ш через 12 констант форм колебаний. 
Для симметричных колебаний (А=1, 2) 
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для дважды вырожденных колебаний (А=3, 4; 5=1, 2) 
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Для определения последних была использована 
система из 12 уравнений, получаемых из 6-ти усло- 
вий ортогональности 
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2-х первых и 2-х вторых условий Эккарта для дваж- 
ды вырожденных колебаний 
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Левые части этих уравнений представляют со- 
бой линейные комбинации от констант форм коле- 
баний, которые необходимо определить, и шести 
квадратичных силовых постоянных, последние бы- 
ли взяты из [4], а правые части равны нулю. Итак, 
система из 12 уравнений была получена с исполь- 
зованием программ, разработанных на языке ана- 
литического программирования МАРЬЕ. Для ре- 
шения этой нелинейной системы в численном ви- 
де была создана программа на языке программиро- 
вания РОКТКАИ, позволяющая решить данную 
систему уравнений методом численного диффе- 
ренцирования. Полученные в результате констан- 
ты форм колебаний приведены в табл. 1. 

Таблица 1. Ненулевые константы форм колебаний молекулы 
АвН з 



Значение 


Значение 


Значение 

/щ. 

Значение 

1]х\ 

0,49427 

кті 

-0,48459 

4узі 

0,29956 

4г41 

0,16554 

/|х2 

-0,29837 

/2x31 

0,16067 

Ьу4\ 

-0,40033 

/4x31 

-0,11608 

4і 

-0,29252 

4x41 

0,48077 

4г31 

0,22123 

4x41 

-0,08686 


Найденные на данном этапе константы форм 
колебаний позволяют восстановить гармонические 
частоты с хорошей точностью. Следовательно, все 
параметры, как г е , а‘, /%, Р ір 1 Ш можно использо- 
вать в качестве нулевого приближения для даль- 
нейших расчетов при определении потенциальной 
функции арсина. Результаты расчета и гармониче- 
ские частоты из [4] приведены в табл. 2. 

Таблица 2. Гармонические частоты молекулы А5Н 3 


Ю А 

Наши результаты, см - ’ 

Результат из [4], см -1 

СО ] 

2189,98 

2190 

®2 

924,01 

924 

Юз 

2201,94 

2202 

0) 4 

1014,97 

1015 


Левые части этих 10-ти уравнений представля- 
ют собой функции от констант форм колебаний, 
все коэффициенты которых определены в числен- 
ном виде, а правые части равны 0 или 1. Два недо- 
стающих уравнения были взяты из перекрестных 
условий на вторые производные от потенциальной 
функции 


Полученные результаты использовались далее 
для определения квадратичной части потенциаль- 
ной функции молекулы путем решения обратной за- 
дачи. При этом была взята потенциальная функция в 
симметрично адаптированных координатах [11]: 

ѵ = к+\Ті р и 5 < 8 і +дЕ А . $ ъ ь X ц» 8 > ^ 5 , , 

^ і) О і ік і ікі 

где Р ір Р і]к , Р т - силовые постоянные; А, 8 Р 8 Ь 8, - 
симметричные координаты. 

Для молекулы арсина симметричные координа- 
ты связаны с естественными, описывающими из- 
менения длин валентных связей и углов между ни- 
ми следующим образом: 
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где а, - изменение угла между /• и связью, г АІ - из- 
менение /-ой связи (/,/, к= 1 , 2 , 3). 

Квадратичная часть потенциальной функции в 
симметрических координатах молекулы типа ХѴ 3 
симметрии С 3ѵ , имеет вид: 

ѵ (2) = \[р п % +ад 2 + 2 /ъад + ^зз + 4)]+ 
+\[гЖ. + й)+ г/^дад + ад )]. 


Следует отметить, что использование симметри- 
ческих координат в теоретических формулах оказы- 
вается не совсем удобным, поэтому важно уметь 
определять связи между параметрами потенциаль- 
ной функции в симметрических и в нормальных ко- 
ординатах, поскольку использование нормальных 
координат более эффективно при теоретических 
расчетах. Процедура перехода от потенциальной 
функции в симметрических координатах к потен- 
циальной функции в нормальных координатах при- 
ведена, например, в [12]. Необходимо отметить, что 
реализация этого перехода для конкретной молеку- 
лы, несмотря на кажущуюся простоту формул, пред- 
ставляет собой довольно сложную задачу, поскольку 
дальнейшие вычисления имеют очень громоздкую 
форму, и проделать их вручную практически невоз- 
можно. Поэтому нами были разработаны програм- 
мы на языке аналитического программирования 
МАРЬЕ для вычислений на компьютере необходи- 
мых аналитических формул. 

В результате проведенных преобразований и 
расчетов мы получили квадратичную часть потен- 
циальной функции молекулы в виде линейной ком- 
бинации шести квадратичных силовых постоянных 
Р п , Е 2Ъ Р а , Р 22 , Р А 4 , Р ЗА , и определили зависимость 
гармонических частот молекулы от Р 9 параметров. 
Коэффициентами при Р і2 являются весьма сложные 
функции от констант форм колебаний молекулы, 
масс атомов, равновесных координат ядер. Следует 
заметить, что все эти величины найдены (табл. 1 ). 
Иными словами, все коэффициенты при Р ѵ могут 
быть определены в численном виде: 

й)!=(0, 170- КЬЕц+0, 164- 1 0 4 Е 22 +0, 106- 10 6 Е 12 ) 1/2 , 

й) 2 =(0,453-10 4 і 7 и +0Л47-10 7 і 7 22 -0,163'10 6 і 7 12 ) 1/2 , 

й) 3 =(0,78 3- 10 2 Е 44 -0,23 1 • 10 5 Е 34 +0, 17 1 • 10 7 Е 33 ) 1/2 , (5) 


й) 4 =(0,154-10 7 і 7 44+0,559-10 5 і 7 3 4+0,508-10 3 7 7 33 ) 1/2 . 

Из анализа уравнений (5), следует, что наиболь- 
ший вклад в частоты со { и щ, описывающие валент- 
ные колебания в молекуле АкН, (как можно было 
ожидать), дают силовые параметры Р п и Е 33 , а в ча- 
стоты <у 2 и ы А — Р 21 и Р ц , поскольку эти колебания 
являются деформационными. 


На первом этапе определения потенциальной 
функции молекулы арсина мы использовали в каче- 
стве начального приближения параметры гармони- 
ческой части Р р полученные в работе [4]. В качестве 
экспериментальных уровней, использовавшихся в 
итерационной процедуре для определения параме- 
тров потенциальной функции, ангармонических х А{ , 
и §„■ параметров, мы использовали информацию о 
26 центрах полос, исследованных в работах [ 6 - 8 ] . В 
том числе 4 фундаментальные полосы ѵ, (А { ), ѵ 2 (4), 
ѵ 3 ( Е) и ѵ А (Е), все первые обертоны 2ѵ, (А 3 ), 2ѵ 2 (Д), 
2ѵ 3 (Д, Е) и 2 ѵ 4 ( Е) (6 полос), дважды возбужденные 
комбинационные полосы ѵ 2 +ѵ 4 (Е) и ѵДѵ 3 ( Е ). В 
этом случае полученные с экспериментальной точ- 
ностью (которая достигает 0,001. ..0,003 см -1 для раз- 
личных полос и различных спектральных диапазо- 
нах) численные значения центров полос позволяют 
нам корректно определить значения гармониче- 
ских частот со А , ангармонических постоянных х А; , и 
наиболее важные параметры резонансных взаимо- 
действий. Обратная задача решалась с использова- 
нием колебательной матрицы гамильтониана, учи- 
тывающей колебательные резонансные взаимодей- 
ствия различных типов: 

Я,ЙГ =ЕИ<%»" 


Здесь суммирование ведётся по всем исследуе- 
мым колебательным состояниям іэ, ѵ. Диагональ- 
ные элементы матрицы имеют вид: 



где со х - значение гармонической частоты; х Х)1 и х Х)К 
- ангармонические постоянные; §„■ - параметры, 
отвечающие за расщепление колебательных термов 
с одинаковыми значениями квантового числа ѵ х , 
Ор г>„; / А - колебательные квантовые числа; (1 = I 
для невырожденных колебательных состояний и 
Д= 2 для дважды вырожденных колебательных со- 
стояний; / А =-и А , -ѵ х +2,...,ѵ х —2,ѵ х . 

Что касается недиагональных матричных эл- 
ементов, то их вид зависит от типа резонанса, кото- 
рый возникает между соответствующими колеба- 
тельными уровнями. При решении нашей задачи 
важно учитывать два типа резонансных взаимодей- 
ствий. Это резонанс Дарлинга-Деннисона, кото- 
рый возникает между валентными колебаниями, 
если Аіэ=+2 и Д/=0, в этом случае недиагональные 
матричные элементы имеют следующий вид: 

Ко = -^Ч(4 + 1)(Ч + 2)(г>з - / 3 2 )] 1/2 , 

если|и)=(и 1 ,...,ц 3 ,...,/з,...),и|и)=(и 1 +2,...,ц 3 -2,...,/ 3 ,...),т.е. 

До, = -2, Аѵ } = +2, А/, = 0 (6) 

и Ко ~ - [°і (°і - 1)((°з + 2 ) 2 - / 2 )] 1/2 , 
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если|и)=(и 1 ,...,Ц 5 ,...,/з,...),и| 0 )=(и 1 - 2 ,...,и 3 + 2 ,...,/з,...),т.е. 

Ац = +2, Аи 3 = -2, А / 3 = 0. (7) 

Здесь параметр квартичной части потенциальной 
функции к { ззз определяет главный вклад резонансно- 
го взаимодействия типа Дарлинга-Деннисона. 

В случае резонанса, который возникает между 
валентными колебаниями, если Дг>=+1 и Д/=±3, 
матричные элементы имеют следующий вид: 

Кв = ± |*іззз [Ц ( и з + к )(°з ± к + 2 )(^з ± к + 4 )] ‘ /2 > 

если |и)=(ц-1,...,ц^+1,...,/з±3,...), и |и)=(ц,..., и,,.. т.е. 

До, =+1, До 3 = -1, А/ 3 = ±3. (8) 

и 

А ии - = + ^ ^ 1333 [(и, + 1)(и 3 + / 3 )(о 3 ± / 3 + 2)(и 3 + / 3 - 2)] 1/2 , 

если |г>>=(г> 1 +1,...,и,-1,...,/ 3 +3,...), и т.е. 

До, =-1, Аи 3 =+1, Д/ 3 = ±3. (9) 

Здесь параметр квартичной части потенциаль- 
ной функции Агцзз определяет главный вклад резо- 
нансного взаимодействия типа Ферми. 


Таблица 3. /у параметры, полученные из экспериментальных 
данных 


Параметр 

Наш результат 

Работа [4] 

А, аДж-А~ 2 

2,83771(279) 

2,810 

Т\ 2 , аДж-А~’ 

0,138* 

0,138 

Тц, аДж 

0,58465(315) 

0,588 

Дз, аДж-А~ ! 

2,88249(190) 

2,839 

Р 34 , аДЖ'А~’ 

-0,047* 

-0,047 

Ян, аДж 

0,67200(257) 

0,669 

Число экспериментальных центров полос 

Число варьируемых параметров 
Среднеквадратичное отклонение, гтз 

26 

14 

0,58 см~’ 


*значение параметра фиксировано значению, взятого из [4] 


Таблица 4. Колебательные спектроскопические параметры 
молекулы АзН 3 , см~' 


Пара- 

метр 

Наш резуль- 
тат 

Резуль- 
тат из [4] 

Пара- 

метр 

Наш резуль- 
тат 

Резуль- 
тат из [4] 

й)і 

2199,29572* 

2190 

*24 

-1,862(987) 

0,03 

С 02 

921,83883* 

924 

Хзз 

-20,033(182) 

-16,38 

СОз 

2218,45970* 

2202 

Х34 

-9,46** 

-9,46 

СО 4 

1017,04996* 

1015 

*44 

-2,106(612) 

-2,42 

*11 

-13,205(235) 

-10,88 

йз 

6,290(155) 

5,39 

Х,2 

-2,06** 

-2,06 

Й4 

-0,36** 

-0,36 

Хіз 

-51,830(330) 

-43,28 

#44 

1,629(204) 

2,12 

*14 

-5,49** 

-5,49 

- 

- 


Хц 

-3,677(645) 

-3,90 

^1133 

49,689(399) 


*23 

-4,84** 

-4,84 

^1333 

-11,9867(481) 



*значения не варьировались, рассчитаны с Е {; параметрами 
из колонки 2 табл. 3, 

**значение параметра фиксировано значению, взятого из [4] 


Таблица 5. Значения центров полос молекулы АзН 3 , см ' 


Полоса 

Из [4] 

Наш результат 

Эксперимент 

Ссылка 

ѵ 2 , А 

910,797 

906,752 

906,752 

[8] 

2 ѵ 2 ,А 

1813,79489 

1806,14959 

1806,14959 

[8] 

2 ѵ 4 , Д 

1985,7889 

1990,9811 

1990,9811 

[7] 

ѵі, А 

2117,824 

2115,260 

2115,164 

[8] 

2ѵз, А 

4205,478 

4167,699 

4166,772 

[8] 

2ѵ*і, А 

4213,888 

4238,046 

4237,701 

[б] 

Ѵ!+2Ѵз, А 

6177,717 

6136,575 

6136,316 

[б] 

Зѵі, А 

6288,19 

6276,00 

6275,83 

[б] 

Зѵ 3 , а 

6307,59 

6365,73 

6365,96 

[б] 

4ѵз, А 

8027,057 

8027,541 

8028,977 

[б] 

4ѵі, А 

8246,246 

8249,19 

8249,52 

[б] 

ѵ 4 , Е 

997,434 

999,225 

999,225 

[7] 

Ѵг+Ѵц, Е 

1908,26191 

1904,11458 

1904,11458 

[8] 

2ѵ 4 , Е 

1994,2690 

1997,4957 

1997,4957 

[7] 

Ѵз, Е 

2124,509 

2126,970 

2126,423 

[8] 

Ѵі+Ѵз, Е 

4199,053 

4168,188 

4167,935 

[8] 

2ѵ 3 , Е 

4227,04 

4248,64 

4247,53 

[б] 

ѵі+2ѵз, Е 

6154,94 

6136,57 

6136,33 

[б] 

2ѵ'і+Ѵз, Е 

6258,30 

6282,70 

6282,36 

[б] 

Зѵз, Е 

6264,47 

6294,67 

6294,71 

[б] 

Ѵі+Зѵз, Е 

8047,712 

8027,547 

8028,969 

[б] 

2ѵі+2ѵз, Е 

8252,45 

8258,13 

8258,38 

[б] 

5 ѵі, А 

9868,4 

9841,0 

9841,4 

[б] 

бѵі, А 

11272,0 

11577,1 

11576,3 

[б] 

4\\+ѵ 3 , Е 

9851,4 

9841,0 

9841,4 

[б] 

5 Ѵ 1 +Ѵ 3 , Е 

11341,5 

11577,1 

11576,3 

[б] 


Квантовые числа г>, и /„ не упомянутые в форму- 
лах (6-9), имеют одно и то же значение в волновых 
функциях |г>) и |г>). 

В результате мы получили ряд из 15 параметров, 
которые представлены в колонке 2 табл. 3 и колон- 
ках 2 и 5 табл. 4 с Іа- доверительным интервалом. 
Полученные параметры воспроизводят экспери- 
ментальные значения центров полос, используе- 
мых в фитинге с Д„ и =0,58 см -1 . В колонках 3 и 6 
табл. 4 приведены для сравнения соответствующие 
параметры, рассчитанные ранее с помощью аЬ іпіГ- 
іо методов [4]. Для того, чтобы проиллюстрировать 
корректность полученных результатов, в колонке 3 
табл. 5 даны значения центров полос, рассчитан- 
ных с нашими параметрами, взятыми из колонок 2 
и 5 табл. 4. Анализ этих данных показывает хоро- 
шее совпадение между экспериментальными и рас- 
считанными значениями центров полос. В колонке 
2 табл. 5 представлены для сравнения значения тех 
же центров полос, рассчитанных с использованием 
параметров, взятых из [4]. 

Следует заметить, что фитинг параметров по- 
тенциальной функции гармонической части осу- 
ществлялся непосредственно к эксперименталь- 
ным значениям центров полос. Таким образом, в 
ходе решения задачи мы нашли, что четыре пара- 
метра Р п , Р 1Ъ Р 3 з и Т 44 определяются однозначно, в 
то время как параметры Р а и Р ц имеют доверитель- 
ные интервалы больше значений самих параме- 
тров. На данном этапе решения задачи параметры 
Р а и Р ц были фиксированы значениям из работы 
[4]. Определение параметров и обуславливает- 


64 







Технические науки 


ся наличием экспериментальной информации. В 
нашем случае, имеющиеся данные о 26 колебатель- 
ных центрах полос позволили нам определить 6 из 
10 параметров х ^ и 2 из 3 параметров § ІГ (табл. 4). 
Аналогично, поскольку часть параметров х Яр не 
удалось определить в силу дефицита эксперимен- 
тальных данных, то их значения были фиксирова- 
ны значениям параметров из работы [4], получен- 
ных из аЬ іпіііо расчетов. Учет резонансных взаимо- 
действий позволил также определить два параме- 
тра квартичной части к ш и к ш . 
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В заключение отметим, что, несмотря на про- 
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ты могут использоваться в дальнейшем как хоро- 
шая основа для коррекции параметров потенциаль- 
ной функции молекулы арсина, полученных с по- 
мощью аЬ іпіііо расчетов. Кроме того, полученная 
информация полезна для дальнейшего спектроско- 
пического анализа слабых высоковозбужденных 
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